Constraints on the high-density nuclear equation of state from the phenomenology of 

compact stars and heavy-ion collisions 
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A new scheme for testing nuclear matter equations of state (EsoS) at high densities using con- 
straints from neutron star (NS) phenomenology and a flow data analysis of heavy-ion collisions 
is suggested. An acceptable EoS shall not allow the direct Urea process to occur in NSs with 
masses below 1.5 Mq, and also shall not contradict flow and kaon production data of heavy-ion 
collisions. Compact star constraints include the mass measurements of 2.1 ± 0.2 Mq (la level) for 
PSR J0751+1807 and of 2.0 ±0.1 M Q from the innermost stable circular orbit for 4U 1636-536, the 
baryon mass - gravitational mass relationships from Pulsar B in J0737-3039 and the mass-radius 
relationships from quasiperiodic brightness oscillations in 4U 0614+09 and from the thermal emis- 
sion of RX J1856-3754. This scheme is applied to a set of relativistic EsoS constrained otherwise 
from nuclear matter saturation properties with the result that no EoS can satisfy all constraints, 
but those with density-dependent masses and coupling constants appear most promising. 

PACS numbers: 04.40.Dg, 12.38.Mh, 26.60.+C, 97.60.Jd 
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I. INTRODUCTION 

The investigation of constraints for the high-density 
behavior of nuclear matter (NM) has recently received 
new impetus when the plans to construct a new acceler- 
ator facility (FAIR) at GSI Darmstadt were published. 
Among others a dedicated experiment for the investiga- 
tion of the phase transition from hadronic matter to the 
quark-gluon plasma (QGP) in compressed baryon matter 
(CBM) shall be hosted, which will study the phenomena 
of chiral symmetry restoration and quark (gluon) decon- 
finement accompanying the transition to the QGP. A firm 
theoretical prediction for the critical baryon densities and 
temperatures of this transition in the QCD phase dia- 
gram as well as the existence and the position of a critical 
point depends sensitively on both the properties of NM 
at high densities and the model descriptions of quark- 
gluon matter in the nonperturbative regime close to the 
hadronization transition. 
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In the present work we apply recently discovered as- 
trophysical bounds on the high-density behavior of NM 
in (3- equilibrium, i.e. neutron star matter (NSM), from 
compact star cooling phenomenology and neutron star 
mass measurements together with information about the 
elliptical flow in heavy-ion collisions (HICs) in order to 
suggest a scheme for testing NM models. This new 
test scheme will be applied to candidates for the NM 
EoS which describe properties at the saturation density 
n s w 0.14 — 0.18 fm~ 3 such as the binding energy per 
particle in symmetric nuclear matter (SNM) a v , the com- 
pressibility K and the asymmetry energy J and charac- 
teristics of large nuclei, such as the neutron skin, surface 
thickness and spin-orbit splitting probing the domain of 
subsaturation densities. In this paper we do not discuss 
the possibilities of various phase transitions, like hyper- 
onization, pion and kaon condensations, quark matter 
etc. 0, Corresponding comments on how their 

inclusion could affect our results are added at the appro- 
priate places. 

While there are several NM models giving a rather 
similar description of the saturation and subsaturation 
behavior they differ considerably in their extrapolations 
to densities above ~2 n Sl the regime which is relevant 
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for NS physics and heavy-ion collisions. Recent progress 
in astrophysical observations and new insights into the 
compact star cooling phenomenology allow us to suggest 
in this paper a test scheme for the high density EoS which 
consists of five elements. 

The first one demands that any reliable nuclear EoS 
should be able to reproduce the recently reported high 
pulsar mass of 2.1 ± 0.2 M Q for PSR J0751+1807, a mil- 
lisecond pulsar in a binary system with a helium white 
dwarf secondary £|. Extending this value even to 2a 
confidence level (j^ ' 5 Mq) means that masses of at least 
1.6 Mq have to be allowed. Thus the EoS should be 
rather stiff to satisfy this constraint. 

The second constraint has recently been suggested in 
Ref. p| and concerns pulsar B in the double pulsar system 
J0737-3039 which has the lowest reliably measured mass 
for any NS to date, namely M = 1.249±0.001 M & @. If 
this star originates from the collapse of an ONeMg white 
dwarf [5J and the loss of matter during the formation of 
the NS is negligible, the baryon number, or equivalently 
the corresponding free baryon mass for the NS, has been 
determined to 1.366 Af© < M N < 1.375 Mq. It turns 
out that this constraint requires a rather strong binding 
of the compact star. A possible baryon loss of up to 1% 
of Mq during the formation of the compact star broadens 
the corresponding baryon mass interval to 1.356 Mq < 
M N < 1.375 Mq. 

The next constraint emerges from recent results of NS 
cooling calculations and population synthesis models 
for young, nearby NSs . Following the arguments in 
Refs. 0,1|, direct Urea (DU) processes, e.g., the neu- 
tron /3-decay n — > p + e~ + v e , produce neutrinos very 
efficiently. The neutrino emissivities for these processes 
even with inclusion of nuclcon superfluidity effects are 
large enough that their occurrence would lead to an un- 
acceptably fast cooling of NSs in disagreement with mod- 
ern observational soft X-ray data in the temperature - age 
diagram. According to these recent analyses, the DU pro- 
cess shall not occur in typical NSs which have masses in 
the range of M typ ~ 1.0 -j- 1.5 Mq, obtained from popu- 
lation syntheses scenarios, see and references therein. 
This constrains the density dependence of the nuclear 
asymmetry energy which should not be too strong. 

The fourth constraint defines an upper bound in the 
mass-radius plane for NSs, derived from quasiperiodic os- 
cillations (QPOs) at high frequencies of the low-mass X- 
ray binary (LMXB) 4U 0614+09 0. For some LMXBs 
there is evidence for the innermost stable circular orbit, 
which if confirmed suggests that the masses of the NSs 
in many of these systems is between 1.8 Mq and 2.1 Mq 
00. 

The fifth constraint comes from a recent analysis of the 
thermal radiation of the isolated pulsar RX J1856 which 
determines a lower bound for its mass-radius relation that 
implies a rather stiff EoS [l3j . 

Finally, we include into the scheme constraints that 
are derived from analyses of elliptic flow data and from 
kaon production in heavy ion collisions. Nuclear colli- 



sions have been described within a kinetic theory ap- 
proach and the results have been compared to experimen- 
tal data for the nucleon flow for densities up to 4.5 x n s 
|14|. From this a region in the pressure-density diagram 
for SNM has been given which defines upper (UB) and 
lower (LB) bounds to the high density EoS and which 
is in accordance with measurements of the elliptic flow. 
Even though 4.5 x n s is below typical central densities 
that correspond to maximum masses of NS configura- 
tions we use the fact that the existence of such a region 
rules out rather stiff and very soft EsoS. 

The outline of this work is the following. In section 
|Ii|we describe a set of modern relativistic nuclear EsoS 
obtained within different approaches. In the section ITTT1 
the test scheme sketched above will be discussed in de- 
tail. This includes the astrophysical constraints from the 
determination of (maximum) NS masses in IIII A II the 
new mass-baryon number test in IIII A2l constraints for 
DU-cooling in IIII A 31 and for the mass-radius relations 
of LMXBs in IIII A 41 as well as the mass-radius relation 
from thermal emission of the isolated NS RX J 1856 in 
IIII A 51 The EoS for SNM at supernuclear densities is 
constrained by HIC experiments from flow data analysis 
in IIII B ll and kaon production in IIII B 21 In Section llVl 
we derive two immediate consequences of this scheme: a 
conjecture about a universal symmetry energy contribu- 
tion to the EoS in /^-equilibrium and a sharpening of the 
flow constraint from HICs using new information about 
the masses of compact stars. A summary of the results 
of this work is given in section together with the con- 
clusions to be drawn from them. 



II. HADRONIC EOS 
A. Model independent description 

There are numerous comparative studies of NM ap- 
proaches for HIC and NS physics applications in which a 
representation of the NM EoS has been employed which 
is based on the nucleonic part of the binding energy per 
particle given in the form 

E(n,f3)=E (n)+f3 2 E s (n), (1) 

where [3 — 1 — 2x is the asymmetry parameter depending 
on the proton fraction x = n p /n with the total baryon 
density n = n n + n p . In Eq. Q the function Eq(ji) is the 
binding energy in SNM, and Eg(n) is the (a)symmetry 
energy, i.e. the energy difference between pure neutron 
matter and SNM. Both contributions Eo(n) and Eg(n) 
are easily extracted from a given EoS for the cases (3 = 
and P = 1, respectively. The parabolic interpolation has 
been widely used in the literature, see e.g. [15j. It proves 
to be an excellent parameterization of the asymmetry 
dependence for the purpose of the present study and we 
will not go beyond it here. Nevertheless, it should be 
mentioned in this context that an exact reproduction of 
a given EoS might require higher order terms than (3 2 
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which have been neglected here. From Eq. (JTJ all zero 
temperature EsoS of NM can be derived by applying sim- 
ple thermodynamic identities 0|. In particular, we ob- 
tain 

e B (n,/3) = nE(n,p), (2) 

P B (n,/3) - n 2 ^-E(n,(3), (3) 
on 

(p 2 T2p-p 2 n-^E s (n) (4) 

for the baryonic energy density e(n) and pressure P(n) 
as well as the chemical potentials of neutron /Lt n (upper 
sign) and proton p p (lower sign), respectively. 

NSM has to fulfill the two essential conditions of j3- 
equilibrium 

Mn = Mp + Me = Mp + Mm i ( 5 ) 

and charge neutrality 

n p - n e - = , (6) 

where // e and /z^ are the electron and muon chemical po- 
tentials, conjugate to the corresponding densities n e and 
riu- In this paper we do not consider phase transitions to 
a deconfincd phase at n > n s . If a first order phase tran- 
sition were allowed a mixed phase could arise in some 
density interval, see 0|- I* 1 general, the local charge 
neutrality condition could be replaced by the global one. 
However, due to the charge screening this density interval 
is essentially narrowed [la . ITflj . The effect of the mixed 
phase on the EoS is also minor. 

Due to Eq. (JSJ the chemical potentials for muons and 
electrons are equal, ji^ = p e so that muons appear in the 
system, once their chemical potential exceeds their mass. 
The EoS for NSM is considered as an ideal mixture of a 
baryonic and a leptonic part, 

e(n,l3) = e B (n,/3) +e e (n,(3) +£ / ,(n,/3) , (7) 
P{n,0) = P B (n,0)+P e (n,0)+P, i {n,0) . (8) 

Under NS conditions one parameter is sufficient for a 
complete description, e.g. the baryochemical potential /x& 
which is conjugate to the conserved baryonic charge. In 
/3-equilibrated NSM and in SNM it is simply equivalent 
to the neutron chemical potential, pi, = p n . Applying 
Eq. Q and Eq. JSJ shows that the electron and muon 
chemical potential can be written as an explicit function 
of baryon density and asymmetry parameter, 

p e (n,0)=A0E s (n). (9) 

Both electrons and muons are described as a massive, 
relativistic ideal Fermi gas. 

With the above relations only one degree of freedom, 
namely the baryon density, remains in charge neutral and 



/3-equilibrated NSM at zero temperature. Within this 
comfortable description actual properties of NM depend 
on the behavior of Eo(n) and Es{n) only. Both can be 
deduced easily from any EoS introduced in the following 
section. 



B. Equations of state applied in this paper 

A wide range of densities up to and above ten times the 
saturation density of NM is explored in the description 
of NSs and HICs. It is obvious that relativistic effects 
are important under these conditions. Consequently, we 
study only nuclear EsoS that originate from relativistic 
descriptions of NM. There are a number of different ap- 
proaches. 

Phcnomcnological models are based on a relativistic 
mean-field (RMF) description of NM with nucleons and 
mesons as degrees of freedom [gjj, El |U The 
mesons couple minimally to the nucleons. The coupling 
strengths are adjusted to properties of NM or atomic 
nuclei. A scalar meson (a) and a vector meson (uj) 
are treated as classical fields generating scalar and vec- 
tor interactions. The isovector contribution is gener- 
ally represented by a vector meson p. In order to im- 
prove the description of experimental data, a medium 
dependence of the effective interaction has to be in- 
corporated into the model. In many applications of 
the RMF model, non-linear (NL) self-interactions of 
the a meson were introduced with considerable success 
IllEllllllllIillillSlllllllll. This approach 
was later extended to other meson fields As an al- 
ternative, RMF models with density-dependent nucleon- 
meson couplings were developed jH HI |M HI EHH . 
They allow for a more flexible description of the medium 
dependence and several parameterizations were intro- 
duced recently. In our study we choose two versions of 
the NL models with self-couplings of the a meson field 
that were used in the simulation of HICs j3^- In the 
parameter set NLp the isovector part of the interaction 
is described, as usual, only by a p meson. The set NLpS 
also includes a scalar isovector meson S that is usually 
neglected in RMF models jljj. It leads to an increased 
stiffness of the neutron matter EoS and the symmetry en- 
ergy at high densities. These particular NL models were 
mainly constructed to explore qualitatively this scalar- 
isovector contribution in the symmetry energy. However, 
they have no non-linearity or density dependence in the 
isovector sector and lead to very, perhaps too, stiff sym- 
metry energies at high densities. The density dependent 
RMF models are also represented here by two parame- 
ter sets |4j. They are obtained from a fit to properties 
of finite nuclei (binding energies, radii, surface thickc- 
nesses, neutron skins and spin-orbit splittings). The pa- 
rameterization DD is the standard approach with con- 
strained rational functions for the density dependence of 
the isoscalar meson couplings and an exponential func- 
tion for the p meson coupling (3{|. In the D 3 C model 
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a 


555 


11.024 


1.4867 


0.19560 


0.42817 0.88233 




783 


13.575 


1.5449 


0.18381 


0.43969 0.87070 


P 


763 


3.6450 


0.44793 







TABLE I: Parameters of the DD-F model as defined in Ref. 
|4*H with a reference density of n rc f = 0.1469 fm 3 . 



additional couplings of the isoscalar mesons to deriva- 
tives of the nucleon field are introduced that lead to a 
momentum dependence of the nucleon self-energies that 
is absent in conventional RMF model |4lj . 

Finally, we present a new parameterization of the RMF 
model with density-dependent couplings that is fitted 
to properties of finite nuclei (binding energies, charge 
and diffraction radii, surface thicknesses, neutron skin in 
208 Pb, spin-orbit splittings) as in Ref. [41j with an addi- 
tional flow constraint (see below) by fixing the pressure of 
SNM to P = 50 MeV fm" 3 at a density of n = 0.48 fm" 3 . 
The density dependence of the a and u> meson coupling 
functions is written as 



Ti(n) = a. 



l + bi(x + dif 



rj(nref)j 



(10) 



l + c^x + d,) 2 ' 

where for the p meson a simple exponential law 

Y p (n) = T p (n rct ) exp[-a p (x - 1)] (11) 

is assumed. The coupling constants Fj(n re f) have been 
fixed at a reference density n ro f. The density dependent 
couplings are functions of the ratio x = n/n IC { with the 
vector density n. The parameters of this parameteriza- 
tion called DD-F are specified in Table [i] For a more 
detailed description of these type of models see Ref. j^J • 

More microscopic approaches start from a given free 
nucleon-nucleon interaction that is fitted to nucleon- 
nucleon scattering data and deuteron properties. In these 
ab initio calculations based on many-body techniques one 
derives the nuclear energy functional from first principles, 
i.e., treating short-range and many-body correlations ex- 
plicitly. A successful approach to the nuclear many-body 
problem is the Brueckncr hole- line expansion. In the rela- 
tivistic Dirac-Brueckner-Hartree-Fock (DBHF) approach 
|42j| the nucleon inside the medium is dressed by the 
self-energy £ based on a T-matrix. The in-medium T- 
matrix which is obtained from the Bethe-Salpeter equa- 
tion plays the role of an effective two-body interaction 
which contains all short-range and many-body correla- 
tions in the ladder approximation. Solving the Bethe- 
Salpeter equation the Pauli principle is respected and 
intermediate scattering states are projected out of the 
Fermi sea. The summation of the antisymmetrized T- 
matrix interactions with the occupied states inside the 



Fermi sphere yields finally the self-energy in Hartree-Fock 
approximation. This coupled set of equations constitutes 
a self-consistency problem which has to be solved by it- 
eration. It is possible to extract the nucleon self-energies 
from DBHF calculations which can be compared with 
the corresponding quantities in phenomenological RMF 
models, but this is not completely unambiguous as dis- 
cussed in Ref. |44|. Here, we use recent results of (asym- 
metric) NM calculations in the DBHF approach with the 
relativistic Bonn A potential in the subtracted T-matrix 
representation |43 l itlliBl l4t| . 

In order to bridge the gap between fully microscopic 
and more phenomenological descriptions that can be ap- 
plied more easily to various systems, it is often useful to 
adjust the parameters of the latter model to results ex- 
tracted from the former method. As an example of this 
approach, we use a nonlinear RMF model (KVR) with 
couplings and meson masses depending on the a- meson 
field 0. The parameters were adjusted to describe the 
SNM and NSM EoS of the Urbana-Argonne group H3 
at densities below four times the saturation density. Ad- 
ditionally, we study also a slightly modified parameter 
set (KVOR) of this RMF model that allows higher maxi- 
mum NS masses. KVR and KVOR models elaborate the 
fact that not only the nucleon but also the meson masses 
should decrease with increasing NM density. Being moti- 
vated by the Brown- Rho scaling assumption, see |48| , and 
the equivalence theorem between different RMF schemes, 
these models use only one extra parameter compared to 
the standard NL RMF model (NL model). 

The nuclear EsoS of these various models can be char- 
acterized by comparing the parameters in the approxi- 
mation of the binding energy per nucleon 
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K' 
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around saturation as a function of the density deviation 
e = (n — n s ) jn s and the asymmetry [3. In this form the 
EoS is characterized at saturation by the binding energy 
ay, the incompressibility K and its derivative, the skew- 
ness parameter K' and by the symmetry energy J and 
the symmetry energy derivative or symmetry pressure 
L for asymmetric NM. In table [H] these parameters are 
given for the models employed in this study. Addition- 
ally, we give the Dirac effective mass = m — £ (at 
the Fermi momentum) in units of the free nucleon mass 
m depending on the scalar self-energy £ of the nucleon. 

There are significant differences between the models. 
The saturation density n s in phenomenological models 
fitted to describe atomic nuclei (DD, D 3 C, DD-F) is in 
the range 0.147 — 0.15 fm -3 . The models NLp, NLp<5 aim- 
ing at a description of low-energy HIC data use the still 
smaller value n s ~ 0.146 fm~ 3 . In contrast to that, the 
"ab-initio" approach (DBHF) shows a saturation density 
that is considerably larger (0.181 fm -3 ). Approxima- 
tions of the Urbana-Argonne type EoS (KVR, KVOR) 
use 0.16 fm~ 3 . The binding energy per nucleon is very 
similar in all models. The incompressibility K spans a 
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Model 


n s 


ay 


K 


K' 


J 


L 


mo 




[fm- 3 ] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


M 


NLp 


0.1459 


-16.062 


203.3 


576.5 


30.8 


83.1 


0.603 


NLp<5 


0.1459 


-16.062 


203.3 


576.5 


31.0 


92.3 


0.603 


DBHF 


0.1810 


-16.150 


230.0 


507.9 


34.4 


69.4 


0.678 


DD 


0.1487 


-16.021 


240.0 


-134.6 


32.0 


56.0 


0.565 


D 3 C 


0.1510 


-15.981 


232.5 


-716.8 


31.9 


59.3 


0.541 


KVR 


0.1600 


-15.800 


250.0 


528.8 


28.8 


55.8 


0.805 


KVOR 


0.1600 


-16.000 


275.0 


422.8 


32.9 


73.6 


0.800 


DD-F 


0.1469 


-16.024 


223.1 


757.8 


31.6 


56.0 


0.556 



TABLE II: Parameters of NM at saturation for various EsoS 
(see text). 



rather wide range from soft (K w 200 MeV) to rather 
stiff (K w 275 MeV). A major difference is found for the 
derivative K' of the incompressibility that is relevant for 
the densities above saturation. Models with parameters 
that are fitted to properties of finite nuclei (DD, D 3 C) 
lead to a negative value of K' with a rather stiff EoS 
at higher densities. It is well known that the ratio of 
the surface tension to the surface thickness is determined 
by the parameters K and K' |49L IHcj. however, the ex- 
act relation depends on the assumption for the shape of 
the surface. In the microscopic DBHF approach and the 
phenomenological models NLp, NLp<5, KVR, KVOR the 
parameter K' is rather large and, correspondingly, the 
EoS of symmetric matter is softer at high densities. The 
DD-F model constructed here is an exception. In this 
parametrization we wanted to satisfy simultaneously the 
description of finite nuclei and the flow constraint that 
requires a soft EoS at high densities leading to a very 
large K' . Correspondingly, the surface properties are not 
optimally well described by the DD-F model with clear 
systematic trends (radii too small for light nuclei and 
too large for heavy nuclei, too small surface thicknesses 
as compared to experimental data). We also remark that 
the parameters of the nonlinear models NLp and NLp<5 
are not representative for conventional NL models that 
are fitted to properties of finite nuclei, e.g., NL3, for 
which one finds K = 271.5 MeV, K' = -203.0 MeV, 
J = 37.4 MeV, L = 100.9 MeV and m D = 0.596 m 

The symmetry energy J is very similar for all models 
with the exception of a slightly larger value in the DBHF 
calculation. Here one has, however, to keep in mind that 
this value is read off at a correspondingly larger density. 
At n = 0.16 fm" 3 DBHF gives a value of J = 31.5 MeV, 
which is in good agreement with the empirical models 
and also with the variational approach of |42| . 

In contrast, the derivatives L of the symmetry energy 
of the various models are spread over a large range. This 
quantity is closely related to the stiffness of the symmetry 
energy at high densities. In order to describe the experi- 



mental neutron skin thicknesses in atomic nuclei a small 
slope of the neutron matter EoS is required 0, |U • 
Models with L < 60 MeV (DD, D3C, KVR, DD-F) fulfill 
this requirement by introducing an effective density de- 
pendence of the p meson coupling to the nucleon which 
goes beyond conventional NL RMF models. A too small 
value for L on the other hand seems to be in conflict with 
data from isospin diffusion in heavy-ion collisions 54] so 
that recently from a combination of these data the limits 
62 MeV < L < 107 MeV have been suggested, see |E3 
and Refs. therein. Only the models NLp, NLpS. DBHF 
and KVOR satisfy this requirement. However, as our em- 
phasis is on high density constraints of the EoS we will 
not elaborate further on this interesting point here but 
remark that it deserves a proper treatment. 

The Dirac effective mass mjj of the nucleon that ap- 
pears in the relativistic dispersion relation of the nucleons 
also shows a large variation in the comparison. In order 
to describe the spin-orbit splitting in atomic nuclei, a 
small value, typically below or around 0.6 m is required. 
Parameter sets with larger values (KVR, KVOR) might 
have a problem in this respect with the construction of a 
proper spin-orbit potential. Larger values of the effective 
Dirac nucleon mass are motivated by fitting the single 
nucleon spectra in nuclei with a large Landau mass 
m* L ~ 0.9-1.0 m. The works [E3] findm^ ~ 0.74-0.82 m 
from the analysis of neutron scattering off lead nuclei. 
The latter values relate to tod ~ 0.7 — 0.8 m |58|. For 
a recent discussion of the momentum and isospin depen- 
dence of the in- medium nucleon mass, see e.g. Ref. |44j . 




-3 -3 -3 

n [f m " ] n [fm ] n [fm " ] 

FIG. 1: The energy per nucleon in SNM Eo(n) (left panel), 
the symmetry energy Es(n) (middle panel) and the energy 
per nucleon in NSM (/3-equilibrated and charge neutral) for 
the investigated models (right panel). 

The variation in the NM parameters is directly re- 
flected in the behavior of the energy per nucleon in SNM 
Eo(n) and of the symmetry energy Eg(n) at densities 
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above saturation as shown in Fig. ^ The various mod- 
els of this study predict considerably different values for 
E (n) and Es(n) at high densities. Under the condition 
of /3-equilibrium, however, the range of binding energy 
per nucleon E(n, (5) shows a much smaller variation than 
expected from Eo(n) and Es(n). This is shown in the 
right panel of Fig. ^ and discussed further in Sect. IV. 

III. CONSTRAINTS ON THE EOS AT HIGH 
DENSITIES 

In this section we will investigate to what extent the 
different EsoS introduced in Sect. [U] fulfill the various 
constraints. We postpone the discussion of the results of 
these tests to section Ivl after two new consequences from 
our analysis are presented in Sect. IV. 



implicitly the radial distribution of the energy density e 
inside the star. In order to solve this set of differential 
equations, one has to specify the EoS, i.e., the relation 
between P and e for which we take the EsoS introduced 
in the previous Section [H] Wc supplement our EsoS de- 
scribing the NSs interior by an EoS for the crust. For 
that we use a simple BPS model [6(j. Due to uncertain- 
ties with different crust models one may obtain slightly 
different mass-radius relations. 

The stellar radius R is defined by zero pressure at the 
stellar surface, P(R) = 0. The star's cumulative gravi- 
tational mass is given then by M = m(R) and its total 
baryon number is N = N(R). In order to solve the TOV 
equations the radial change of the pressure P starting 
with a given central value at radius r — has to be 
calculated applying, e.g., an adaptive Runge-Kutta algo- 
rithm. 



A. Constraints from compact stars 

1. Maximum mass constraint 

Measurements of "extreme" values, like large masses 
or radii, huge luminosities etc. as provided by compact 
stars offer good opportunities to gain deeper insight into 
the physics of matter under extreme conditions as pro- 
vided by compact stars. Recent measurements on PSR 
J0751+1 807 imply a pulsar mass of 2.1 ± 0.2 M Q 
(first error estimate with la confidence, second in brack- 
ets with 2a confidence) U which is remarkably heavy 
in comparison to common values for binary radio pul- 
sars (M BRP = 1.35 ± 0.04 M This special result 
constrains NS masses to at least 1.6 Mq (2a confidence 
level) or even 1.9 Mq within the la confidence level. 

The mass and structure of spherical, nonrotat- 
ing stars, to which we limit ourselves in this pa- 
per, is calculated by solving the Tolman-Oppenheimer- 
Volkov(TOV)-equation, which reads as 



dP(r) 
dr 



G[e{r) + P(r)][m(r) + 47rr 3 P(r)] 
r[r — 2Gm(r)] 



(13) 



where the gravitational mass m(r) inside a sphere of ra- 
dius r is given by 



m(r) — 4-7T / dr' r' 2 e(r') 



(14) 



which includes the effects of the gravitational binding 
energy. The baryon number enclosed by that sphere is 
given by 



N(r) = An 



dr' r' n(r') 



2Gm(r') 



(15) 



with n(r) being the baryon density profile of the star. 
Eq. i|13|) describes the gradient of the pressure P and 
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FIG. 2: Mass versus central density for compact star con- 
figurations obtained by solving the TOV equations 11311 and 
11411 for all EsoS introduced in Subsect. 2.2. Crosses denote 
the maximum mass configurations, filled dots mark the crit- 
ical mass and central density values where the DU cooling 
process becomes possible. According to the DU constraint, it 
should not occur in "typical NSs" for which masses are ex- 
pected from population synthesis @ to lie in the lower grey 
horizontal band. The dark and light grey horizontal bands 
around 2.1 Mq denote the la and 2a confidence levels, re- 
spectively, for the mass measurement of PSR J0751+1807 0. 

The resulting NS masses as a function of their cen- 
tral density for the different EsoS are given in Fig. [3 to- 
gether with the mass range of typical NSs and the limits 
from PSR J0751+1807. Also shown in this figure are the 
points on the respective curves where the DU process be- 
comes possible, as further discussed in subsection llll A 31 

The maxima of the mass-central density relations are 
easily determined then and summarized in Table IIIII for 
the EsoS investigated in this work. As can be seen none 
of these values falls below the 2a mass limit of 1.6 Mq, 
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whereas the la mass limit of 1.9 Mq would exclude NLp 
and NLp<5, while marginally excluding KVR. Thus the 
ability of this first and rather trivial test to exclude a 
given EoS demands a high accuracy of observations. A 
more stringent test could be achieved with decreasing er- 
ror estimates or the observation of at least one pulsar that 
is still more massive than PSR J0751+1807. We point 
out that if a pulsar with a mass M > 2.1 Mq is observed 
in the future, this will imply serious restrictions on the 
viable EoS, see Fig. Within the set of EsoS tested by 
us, only DD, D 3 C and DBHF would survive. Moreover, 
the maximum mass constraint is closely related to the 
flow constraint. This point will be further investigated 
within subsection IIV Bl 
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TABLE III: Maximum star masses, corresponding central 
densities and the fulfillment of the strong (la) and weak (2a) 
maximum mass constraint, as well as the gravitational mass 
- baryon number constraint for Pulsar B in J0737-3039 5] 
without and with a mass loss of 0.01 Mq. Fulfillment (viola- 
tion) of a constraint is indicated with +( — ) and a marginal 
result is rated with o. 



2. Gravitational mass - baryon number constraint 

Recently, it has been suggested in Q that pulsar B 
in the double pulsar system J0737-3039 may serve to 
test models proposed for the EoS of superdense nuclear 
matter. The system J0737-3039 consists of a 22.7 ms 
pulsars J0737-3039A (pulsar A) H3, and a 2.77 ms pul- 
sar companion J0737-3039B (pulsar B) (6^, orbiting the 
common center of mass in a slightly eccentric orbit of 
2.4 hours duration. One of the interesting characteris- 
tics of this system is that the mass of pulsar B is merely 
1.249 ± 0.001 M Q which is the lowest reliably mea- 
sured mass for any NS to date. Such a low mass could 
be an indication that pulsar B did not form in a type-H 
supernova, triggered by a collapsing iron core, but in a 
type-I supernova of an ONeMg white dwarf [j| driven hy- 
drostatically unstable by electron captures onto Mg and 



Ne. The well-established critical density at which the 
collapse of such stars sets in is 4.5 x 10 9 g/cm 3 corre- 
sponding to an ONeMg core whose critical baryon mass 
is Mm — N u ~ 1.37 Mq, where the atomic mass unit 
u = 931.5 MeV has been used Q to convert the baryon 
number to baryon mass. Assuming that the loss of mat- 
ter during the formation of the NS is negligible, a pre- 
dicted baryon mass for the NS of Mm = 1.366 — 1.375 Mq 
was derived in 5]. This theoretically inferred baryon 
number range together with the star's observed gravi- 
tational mass of M — 1.249 ± 0.001 Mp^ may represent 
a most valuable constraint on the EoS |5j, provided the 
above key assumption for the formation mechanism of 
the pulsar B is correct. Then any viable EoS proposed 
for NSM must predict a baryon number in the range 
1.366 S; Mm S; 1.375 Mq for a NS whose gravitational 
mass is in the range M = 1.249 ± 0.001 Mq. None of the 
EsoS tested in this work satisfies this strong constraint. 
The authors of @ discussed caveats such as baryon loss 
and variations of the critical mass due to carbon flashes 
during the collapse. This constraint requires a very pre- 
cise calculation of the baryon number, e.g. a lowering of 
Mm by 1% changes the outcome of this test significantly. 
Since the simulation of e-capture supernovae and the evo- 
lution of their progenitors is still a work in progress, more 
interesting results are expected in the near future. The 
final value and accuracy of the baryon number of J0737- 
3039 are therefore highly important. The result of such 
calculations is shown in Fig. [3] and summarized in Ta- 
ble [ffi] Finally we point out that this constraint is criti- 
cally based on the assumption of the formation scenario 
for pulsar B. If this turns out to be incorrect the con- 
straint has to be abandoned. 



3. Direct Urea constraint 

The maximum mass constraint seems to have, at least 
for the EsoS investigated in this paper, a rather small 
exclusion potential. The M — Mm criterion, however, 
would provide more stringent limits only if the assumed 
formation mechanism of pulsar B in RX J0737-3039 and 
the neglect of mass loss prove to be valid. Adding the 
DU criterion will be seen to improve this scheme. 

If the proton fraction x = n p /(n p + n n ) exceeds a crit- 
ical value Xdu the DU process n —> p + e~ + v e be- 
comes operative. An estimate of this DU-threshold fol- 
lows from the triangle inequality for momentum conser- 
vation where the moduli of the momenta are given by the 
neutron, proton and electron Fermi momenta pFi- The 
typical neutrino energy of the order of the temperature 
T is small and can be neglected. In quasi-equilibrium 
n — > p + er + v e implies that pF n < Pf + Pf c - From the 



charge neutrality condition n p = n e 
the DU-threshold xdu a s 



1 



xdu = 



one easily finds 



(16) 
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FIG. 3: Relation between gravitational mass M and baryon 
mass Mn (both in units of the solar mass Mq) of NSs for the 
EsoS discussed in this work. The filled rectangle denotes the 
constraint derived for Pulsar B in the double pulsar J0737- 
3039 [3 • The empty rectangle demonstrates the change of 
the constraint when the assumed loss of baryon number in 
the collapse amounts to 1% of the solar value. 



where x e = n e /(n e J rn^) is the leptonic electron fraction. 
Since this depends on the symmetry energy, the DU- 
threshold is model dependent. For x e — 1 (no muons) 
this formula reproduces the muon-free threshold value of 
11.1% 63]. In the limit of massless muons, which is ap- 
plicable for high densities (x e = 1/2) one finds an upper 
limit of Xdu = 14.8%. In Fig.^the proton fraction as a 
function of density is shown for the different models, to- 
gether with the DU-threshold value xdu, given as a band 
for all the models. As can be seen this threshold can be 
reached for a wide range of densities depending on the 
EoS. For some models (DD, D 3 C, DD-F) it does not take 
place at all. The model-dependent DU-threshold occurs 
at a corresponding critical baryon density. Setting this 
as a star's central density results in a DU-critical star 
mass Mdxj. These critical densities and DU-masses are 
marked in Fig. ^ as dots on those model curves, where 
the limit is reached. Every star with a mass only slightly 
above Mdu will be efficiently cooled by DU-processes and 
very quickly becomes almost invisible for thermal detec- 
tion [3 ■ Nucleon superfluidity which suppresses the cool- 
ing rates has been included. Values of the pairing gaps 
used in the literature have been used and then varied to 
check the model dependence of the result |(34[ ■ Table IIVI 
summarizes these DU critical masses for all models. The 
DU constraint is fulfilled by the DD, D 3 C and DD-F EoS 
models which are not affected by the DU process at all 
and by KVR, KVOR which are affected for masses higher 
than the limit for "typical NS" of 1.5 Mq obtained from 
population synthesis models Q . As a weaker constraint 



FIG. 4: Proton fractions x = n p /(n n + n p ) for different 
EsoS. The band labeled with xou frames all threshold curves 
obtained for the investigated models. According to Eq. lib! 
the range of possible threshold values varies between 11.1% 
and 14.8% depending on the muon fraction. 



we also use Muu > 1-35 Mq. This follows from both 
the population synthesis and the mass measurement of 
binary radio pulsars. If the DU process were allowed for 
Mdu < 1-35 Mq it would affect most of the NS pop- 
ulation. It should, however, not be expected that the 
objects observed in X-rays were some exotic family of 
NSs rather than typical NSs. DBHF and both NL mod- 
els do not pass the DU test. They have a DU threshold 
mass Mdu < 1-1 Mq- Note also that only NSs with 
M > 1.1 Mq are produced in the standard scenario of 
NS formation in type II supernova explosions |65j |. 



4- Mass-Radius relation constraint from LMXBs 

The kilohertz quasi-periodic brightness oscillations 
(QPOs) seen from more than 25 NS X-ray binaries con- 
strain candidate high-density EsoS because there are fun- 
damental limits on how high-frequency such oscillations 
can be. A pai r of such QPOs is often seen from these sys- 
tems (see 166] for a general review of properties). In all 
currently viable models for these QPOs, the higher fre- 
quency of the QPOs is close to the orbital frequency at 
some special radius. For such a QPO to last the required 
many cycles (up to ~100 in some sources), the orbit must 
obviously be outside the star. The orbit must also be out- 
side the innermost stable circular orbit (ISCO), because 
according to the predictions of general relativity, inside 
the ISCO gas or particles would spiral rapidly into the 
star, preventing the production of sharp QPOs. This 
implies [ToL l6?| that observation of a source whose maxi- 
mum QPO frequency is f max limits the stellar mass and 
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radius to 



M < 2.2 M Q (1000 Hz/zw)(l + 0.75j) 
R < 19.5 km(1000 Hz/i/ roax )(l + Q.2j) . 



(17) 



Here j = cJ/GM 2 (where J is the stellar angular mo- 
mentum) is the dimensionless spin parameter, which is 
typically 0.1-0.2 for these systems. There is also a limit 
on the radius for any given mass. 

These limits imply that for any given source, the ob- 
served f max means that the mass and radius must fall 
inside an allowed "wedge". Therefore, any allowed EoS 
must have some portion of its corresponding mass-radius 
curve fall inside this wedge. The wedge becomes smaller 
for higher i/ max , therefore the highest frequency ever ob- 
served (1330 Hz, for 4U 0614+091; see HH) places the 
strongest of such constraints on the EoS. Note, though, 
that another NS could in principle have a greater mass 
and thus be outside this wedge, but an EoS ruled out 
by one star is ruled out for all, since all NS have the 
same EoS. As can be seen from Fig. the current con- 
straints from this argument do not rule out any of the 
EoS we consider. However, because higher frequencies 
imply smaller wedges, future observation of a QPO with 
a frequency ~ 1500 — 1600 Hz would rule out the stiffest 
of our EoS. This would therefore be a complementary re- 
striction to those posed by RX J1856. 5-3754 (discussed 
below) and the implied high masses for some specific NSs, 
which both argue against the softest EoS. 

If one has evidence for a particular source that a given 
frequency is actually close to the orbital frequency at 
the ISCO, then the mass is known (modulo slight uncer- 
tainty about the spin parameter). This was first claimed 
for 4U 1820-30 |69|. but complexities in the source phe- 
nomenology have made this controversial. More recently, 
careful analysis of Rossi X-ray Timing Explorer data for 
4U 1636-536 and other sources [Tij has suggested that 
sharp and reproducible changes in QPO properties are 
related to the ISCO. If so, this implies that several NSs 
in low-mass X-ray binaries have gravitational masses be- 
tween 1.9 M Q and possibly 2.1 Mq 11]. In Fig. we 
display the estimated mass 2.0 ± 0.1 Mq for 4U 1636- 
536, which would eliminate NLp and NLpS as the softest 
proposed EoS even in the weak interpretation, and allow 
only DBHF, DD and D 3 C in the strong one, see Tab.HVl 



5. Mass-Radius relation constraint from RX J 1856 

After the discovery of the nearby isolated NS RX 
J1856.5-3754 (hereafter short: RX J1856) the analysis 
of its thermal radiation using the apparent blackbody 
spectrum with a temperature — 57 eV [70| yielded 
a lower limit for the photospheric radius Roo of this ob- 
ject. The distance of RX J1856 was initially estimated 
to be 60 pc. Since R^ crucially depends on this quan- 
tity a very small value of R^ » 8 km was derived which 
could not have been explained even with RX J 1856 be- 
ing a self-bound strange quark star [70j . The true stellar 



radius R is given by R x = R(l — R/Rs)^ 1 ^ 2 , with the 
Schwarzschild radius Rs — 2GM/R. New measurements 
predict a distance of at least 117 pc, which results in 
Roo = 16.8 km and turns RX J1856 from the formerly 
smallest known NS into the largest one [l3| . The result- 
ing lower bound in the mass radius plane is shown in 
Fig- El There are three ways to interpret this result: 

A) RX J1856 belongs to compact stars with typical 
masses M ~ 1.4M Q and would thus have to have 
a radius exceeding 14 km (see Fig- EJ ■ None of the 
examined EsoS can meet this requirement. 

B) RX J1856 has a typical radius of R ~ 12 - 13 
km, implying that the EoS has to be rather stiff 
at high density in order to allow for configurations 
with masses above ~ 2 Mq. In the present work 
this condition would be fulfilled for DBHF, DD and 
D 3 C. This M > 1.6 Mq explanation implies that 
the object is very massive and it is not a typical 
NS since most of NSs have M < 1.5 Mq, as follows 
from population synthesis models. 

C) RX J1856 is an exotic object with a small mass 
~ 0.2 Mq, which would be possible for all EsoS 
considered here. No such object has been observed 
yet, but some mechanisms for their formation and 
properties have been discussed in the literature [7l| . 




FIG. 5: Mass-Radius constraints from thermal radiation of 
the isolated NS RX J1856. 5-3754 (grey hatched region) and 
from QPOs in the LMXBs 4U 0614+09 (green hatched area) 
and 4U 1636-536 (orange hatched region) which shall be re- 
garded as separate conditions to the EsoS. For the mass of 
4U 1636-536 a mass of 2.0 ± 0.1 Mq is obtained, so that the 
weak QPO constraint would exclude the NLp and NLp<5 EsoS 
whereas the strong one would leave only DBHF, DD and D 3 C. 

It cannot be excluded, however, that the distance mea- 
surement could be revised by a future analysis. If the 
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distance would turn out to be smaller than the present 
value, then this constraint would have no discriminative 
power any more since all EsoS could possibly fulfill it. 
Should a revised distance value be larger than the present 
one, then only the exotic low-mass star interpretation 
would remain which again is possible for all (not self- 
bound) EsoS but which would raise the question about 
the formation scenario for such a diffuse low-mass ob- 
ject. Certainly this explanation of the puzzling object 
would no longer qualify RX J1856 as an object to test 
the high-density nuclear EoS. 
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TABLE IV: Critical compact star mass for the occurrence of 
the DU cooling process with the corresponding central den- 
sity, the criterion of the DU constraint and the M(R) con- 
straints from the isolated NS RX J1856. 5-3754 and the low- 
mass X-ray binary 4U 1636-536 with its upper (u) and lower 
(1) mass limits, see text. Symbols are defined in Tabic UTTl 



Finally we want to emphasize another problem. Com- 
paring Fig. 13 and Fig. [5] we observe that models produc- 
ing a smaller radius (KVR, DD-F, DBHF) better accom- 
modate the M — Mm constraint and those having a larger 
radius (D 3 C, DD) fulfill the RX J1856 constraints better. 
Out of the EsoS tested in this work only DBHF could sat- 
isfy these constraints simultaneously and it would be a 
rather challenging task to resolve the problem of this EoS 
with the DU constraint. Further comments are given in 
the discussion below. 



B. Constraints from heavy-ion collisions 

1. The flow constraint 

The flow data analysis of dense SNM probed in HICs 
fl4| reveals a correlation to the stiffness of the EoS which 
can be formulated as a constraint to be fulfilled within 
the testing scheme introduced here. 




FIG. 6: Pressure region consistent with experimental flow 
data in SNM (dark shaded region). The light shaded region 
extrapolates this region to higher densities within an upper 
(UB) and lower border (LB). 



The flow of matter in HICs is directed both forward 
and perpendicular (transverse) to the beam axis. At 
high densities spectator nucleons may shield the transver- 
sal flow into their direction and generate an inhomoge- 
neous density and thus a pressure profile in the transver- 
sal plane. This effect is commonly referred to as elliptic 
flow and depends directly on the given EoS. An analysis 
of these nucleon flow data, which depends essentially only 
on the isospin independent part of the EoS, was carried 
out in a particular model in Ref . [l4| . In particular it was 
determined for which range of parameters of the EoS the 
model is still compatible with the flow data. The region 
thus determined is shown in Fig. as the dark shaded 
region. Ref. then asserts that this region limits the 
range of accessible pressure values at a given density. For 
our purposes we extrapolated this region by an upper 
(UB) and lower (LB) boundary, enclosing the light shade 
region in Fig. H3 

Thus the area of allowed values does not represent ex- 
perimental values itself, but results from transport calcu- 
lations for the motion of nucleons in a collision [l4|. Of 
course, it seems preferable to repeat these calculations for 
each specific EoS, but this would not be a manageable 
testing tool. Therefore we adopt the results of ref. [14( as 
a reasonable estimate of the preferable pressure-density 
domain in SNM. Its upper boundary is expected to be 
stable against temperature variations |72j |. The impor- 
tant fact is that the flow constraint probes essentially 
only the symmetric part of the binding energy function 
E {n). 

Following Ref. Q the constraint arises for a density 
window between 2 and 4.5 times saturation density n s . 
One has, however, to keep in mind that at high densi- 
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ties this constraint is based on flow data from the AGS 
energy regime (£i a b ~ 4 — 11 AGeV). At these ener- 
gies a large amount of the initial bombarding energy is 
converted into new degrees of freedom, i.e., excitations of 
higher lying baryon resonances and mesons, which makes 
conclusions on the nuclear EoS more ambiguous than at 
low energies. Nevertheless, the analysis of 01 provides 
a guideline for the high density regime which we believe 
to be reasonable. 

As can be seen in Fig. this last constraint is well 
fulfilled by the KVR, KVOR, NLp and NLp<S models. 
For the latter two models this is rather obvious since 
they have already been tested to reproduce flow data. 
The constraint is satisfied for densities below 3 n s by 
DBHF. When comparing our models with the flow con- 
straint below, we separate regions of SIS (n < 3n s ) and 
AGS (n > 3n s ) energies considering weak and strong 
flow constraints. DD and D 3 C, which fulfilled the DU 
constraint well, are significantly above the demanded re- 
gion. We want to emphasize that the DD-F model was 
constructed in this paper to pass this test. It is based on 
a reparametrization of the DD model, in order to satisfy 
the introduced test scheme in most points. 



IV. CONSEQUENCES 
A. Universal symmetry energy conjecture 




n [fm ] 



n [fm ' 



2. Constraints from subthreshold kaon production 

K + mesons were suggested as promising tools to probe 
the nuclear EoS, almost 20 years ago [73]. The first 
channel to open in order to produce a K + meson is 
the reaction NN — > NAK + which has a threshold of 
Eub = 1.58 GeV kinetic energy for the incident nucleon. 
When the incident energy per nucleon in a heavy ion re- 
action is below this value one speaks of subthreshold kaon 
production. This process is particularly interesting since 
it ensures that the kaons originate from the high den- 
sity phase of the reaction. The missing energy has to be 
provided either by the Fermi motion of the nucleons or 
by energy accumulating multi-step reactions. Both pro- 
cesses exclude significant distortions from surface effects 
if one goes sufficiently far below threshold. In combi- 
nation with the long mean free path subthreshold K + 
production is an ideal tool to probe compressed NM in 
relativistic HICs, see Ref. [74| for a recent review. 

Within the last decade the KaoS Collaboration at GSI 
has performed systematic measurements of the K + pro- 
duction far below threshold At subthreshold 
measurements which range from 0.6 to 1.5 AGeV labo- 
ratory energy per nucleon compressions of two to max- 
imally three times n s are reached. Transport calcula- 
tions have demonstrated that subthreshold K + produc- 
tion provides a suitable tool to constrain the EoS of SNM 
at supersaturation densities j7J, l?8l l79| . The theoretical 
analysis of the data implies a soft behavior of the EoS 
in the considered density range consistent with the flow 
constraint at moderate densities (n < 3n s ) and supports 
DBHFjNLp, NLpS, KVR, KVOR and DD-F EsoS in 

snm [13, Hil H3l . 



FIG. 7: Density dependence of the asymmetry contribution 
to the energy per particle (left panel) and of the proton frac- 
tion (right panel) in NSM. Encircled curves correspond to 
EsoS that violate the DU-constraint. 

Investigating the onset of DU processes in section 
1111 A 31 has shown that the DU threshold for the inves- 
tigated models can be reached for rather small baryon 
densities slightly below 2n s for NLp, NLpJ, DBHF or, as 
the most extreme opposite, not at all for DD, D 3 C and 
DD-F. Eq. llofl states that the threshold xdu depends 
on the electron-muon ratio. The electron and muon den- 
sities are determined by their chemical potentials. In 
/^-equilibrium /i e = /i M is given in turn by Eq. as a 
function of the asymmetry parameter (3 and the symme- 
try energy E$- The resulting proton fraction x, shown 
on the right hand side of Fig. \7\ mainly maps the topo- 
logical behavior of Es(n), see Fig.^ As a rule of thumb 
therefore a large proton fraction x is attained for stiff 
symmetry energies Es(n). The NLp, NLpS and DBHF 
models confirm this rule well and the symmetry energy 
used by these models can be sorted out for contradicting 
the present cooling phenomenology, as described in sec- 
tion IIIIA3I Fig. H illustrates both DU- and maximum 
mass constraint. 

The asymmetry contribution /3 2 Es(n) to the energy 
per nucleon (left panel of Fig. [TJl only shows a marginal 
dependence for different EsoS when compared to differ- 
ences in the energy per nucleon of SNM. These form 
a narrow band which allows two important statements. 
First, the behavior of (3 2 Es{n) is to good approximation 
universal for all EsoS which pass the DU-constraint. The 
second conclusion regards the influence of the symmetry 
energy on the mass of NS. Here we find that due to the 
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above universal behavior of f3 2 Es(n) the mass of a star 
is also dominated by the behavior of E (n). In other 
words it seems not very likely that it is possible to in- 
fer any essential properties of Es{n) only from NS mass 
observations, even if Eq(ji) would be perfectly known. 
This point emphasizes the importance of a more detailed 
understanding of the cooling behavior of compact stars 
as an effective tool for probing Es(n) at high densities 
beyond saturation. 

B. Sharpening the flow constraint 

As shown in sections IIII A II and IIII B II the flow con- 
straint was more restrictive on the EoS than the maxi- 
mum mass constraint due to a large estimated error for 
the mass of PSR J0751+1807. But the flow constraint 
itself has uncertainties, represented by the region of pos- 
sible pressures in Fig. EJ 

Thus an EoS which is in accordance with the flow con- 
straint might still violate the maximum NS mass con- 
straint. This is demonstrated in Fig. [Sj where in partic- 
ular the lower boundary of the limiting region in Fig. H3 
was extrapolated to construct an artificial EoS (LB) in 
order to obtain the mass-density relation for correspond- 
ing compact star configurations . 
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FIG. 8: Mass versus central density for compact star config- 
urations, calculated using the UB and LB extrapolations of 
the flow constraint boundaries from Fig. |SJ Eo(n), together 
with different symmetry energies Es(n) not violating the DU- 
constraint, see Fig.Q The error bars illustrate the maximum 
deviation resulting from choosing different symmetry energies 
Es(n). The gray horizontal bars denote the expected mass 
of PSR J0751+1807 including for la and 2a confidence inter- 
vals, resp., whereas the vertical bars limit the density region 
covered by the flow constraint. 

The resulting mass curve in Fig. is far from reaching 
even the weak mass constraint of 1.6Mq originating from 



the lower 2a bound on the mass of PSR J0751+1807. 
It also cannot accommodate the well-measured mass 
-^51913+16 = 1.4408±0.0003M Q This figure clearly 
demonstrates that the lower bound (LB) in Fig. does 
not satisfy the maximum mass constraint and should be 
shifted upwards, thus narrowing the band of the flow con- 
straint. 

The maximum mass constraint demands a certain stiff- 
ness of Eo(n) in order to obtain sufficiently large max- 
imum NS masses. The small influence of Es{n) on the 
NS mass can be well recognized on Fig. [SJ too. Two 
different symmetry energies, necessary to describe NSM, 
were taken from the investigated EsoS. They were cho- 
sen in accordance with the DU-constraint and gave the 
largest (DD-F) and smallest (D 3 C) contribution to the 
binding energy at n = 1 fm -3 . The resulting deviations 
of the NS masses are shown as error bars on the curves 
in Fig. EI R results in a maximum difference of less than 
0.2 Mq for the mass curves corresponding to LB. The 
same was done for an artificial EoS extrapolating the up- 
per boundary (UB). Here the largest error estimate of 
approximately 0.1 Mq is even smaller. The maximum 
mass for UB of about 2.0 Mq again fulfills well the cor- 
responding constraint. 

Comparing the flow constraint in Fig. Eland the mass- 
radius constraint in Fig. [5] we see that none of the EsoS 
we tested satisfies both constraints. Only DBHF is able 
to satisfy a weak flow constraint (for n < 3n s ) and the 
RX J1856 constraint under the assumption of an object 
with M > 1.6 Mq. Thus we again emphasize a problem 
for the RX J1856 constraint with the joint fulfillment of 
other constraints as the M—M^ and the flow constraints. 



V. SUMMARY AND CONCLUSIONS 

The task we intended with this work, developing a test 
scheme for the nuclear EoS by the present phenomenol- 
ogy of dense NM in compact stars and heavy-ion colli- 
sions, is satisfactorily completed at this point. Applying 
this scheme to specific EsoS offers some interesting in- 
sights which indicate that astrophysical measurements 
might become more important for the interpretation of 
terrestrial measurements than presently accepted. We 
have summarized the results of all suggested tests per- 
formed on our choice of relativistic, high-density EsoS in 
Table IVI which reveals the discriminative power of their 
combined application in a broad region of densities and 
isospin asymmetry. We want to point out here, however, 
that each model was derived to describe a restricted re- 
gion in the (n, 0)- plane and was not necessarily meant 
to describe a broader region. In Tab. [V]we rate the per- 
formance of the models when applied nevertheless in a 
very wide (n, (3)- region. 

Due to its sensitivity to different contributions to the 
energy this scheme motivates the necessity for changes 
in several EsoS if one wants to apply them to the whole 
available (n, (3) interval although they well describe prop- 
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TABLE V: Summary of results for the suggested scheme of 
tests. Non separated columns show the results for a strict 
(left) and weakened (right) interpretation of the correspond- 
ing constraint. The last column gives the total number of 
fulfilled tests in the suggested scheme. Symbols are defined 
in Table Em 



erties at the saturation density and a specific region of 
the (n, (3) plane. 

In particular Es(n), the contribution of the symme- 
try energy to the total energy, is probed by the DU- 
constraint. It states that the DU process would cool 
NSs much too fast, so that it should occur for stars with 
masses greater than w 1.35 — 1.5M© only. If it would 
affect stars with masses below this limit the DU-process 
would affect most of the known NS population. We have 
shown that only EsoS with a rather soft symmetry energy 
at high density fulfill the DU-constraint. 

The symmetric matter contribution to the energy per 
nucleon Eo(n) should sufficiently describe the elliptic 
flow. Adding the results of 01 to the test scheme, very 
soft EsoS are allowed. The latter, however, can be sorted 
out by the maximum mass constraint. As a result, these 
two combined constraints limit the stiffness of a reliable 
EoS to be rather low (flow), but not too low at n < An s 
and rather high for higher densities (maximum mass). 

We want to stress that out of the tested EsoS only 
DBHF passes simultaneously the gravitational binding 
(M — Mm) as well as both mass-radius (M — R) tests. 
The M — Mm constraint to be fulfilled requires a smaller 
radius of the M ~ 1.25 M Q star, whereas the M — R 
test from RX J1856 favors substantially larger star radii, 
at least for M < 1.4 M©. This contradictory situation 
would be resolved when RX J 1856 is a star with a large 
mass or when an EoS would fulfill the M — Mm constraint 
and nevertheless assign large radii to NS with typical 
masses. Such an EoS would be qualitatively different 
from the ones we investigate here. 

The whole scheme left three of eight model EsoS, 



namely DD-F, KVR and KVOR as most effective within 
a broad (n,/3) region under consideration. The DD-F 
model explicitly fits properties of finite nuclei, especially 
the neutron skin thickness of 208 Pb that implies a small 
value of L. The KVR model yields a similar value of L. 
It might, however, have problems with respect to isospin 
diffusion in heavy ion collisions since this small L does 
not fit the constraint deduced in Ref. |85j and [55j. The 
KVOR model fulfills this latter constraint. Here we point 
out that both KVR and KVOR were not applied to finite 
nuclei. Thus it would be a challenge to apply such mod- 
els to finite nuclei in the future. All the models except 
DD-F demonstrate their predictive power within a broad 
(n, f3) region whereas the DD-F model (as a modified DD 
model) has been constructed in the present work in or- 
der to fulfill the flow constraint in addition to constraints 
from saturation and finite nuclei common to DD models. 
In contrast to the phenomenological RMF models, DBHF 
is an ab initio approach without room left for the read- 
justment of free parameters. But correlations beyond the 
ladder approximation are not taken into account. In a 
certain sense, they are included, although hidden in the 
fitted parameters, in the phenomenological approaches. 
However, an interesting aspect would be to perform cal- 
culations for different types of free space nucleon-nucleon 
interactions. In particular the CD-Bonn potential j8|| 
which accounts more precisely for the isospin dependence 
of the nuclear forces than Bonn A (used here) would 
be appropriate for future investigations. Another point 
would be the explicit inclusion of hyperonic degrees of 
freedom which may have a significant impact on the NS 
matter EoS at high densities (depending on yet badly 
known nucleon- hyperon interaction). This could open a 
possibility for the DBHF and other EsoS to satisfy ap- 
propriately the DU constraint. 

Beside the scheme's good overall selectivity the joint 
application of different constraints might give new in- 
teresting insights. One of these is the universal behav- 
ior of the contribution (3 2 Es{n) to the binding energy 
in NSM we observed for all EsoS that fulfill the DU- 
constraint. Then it seems to us that the flow constraint 
limits the maximum mass of NSs to values around or not 
much about the expected mass of PSR J0751+1807 with 
M = 2.1 Mq, which also coincides with the upper mass 
limit for 4U 1636-536. To verify this suggestion, a more 
detailed analysis, similar to that shown in [87|], has to be 
performed. 

Next we want to emphasize that the maximum mass 
constraint as a result of astrophysical measurements fur- 
ther limits the pressure-density-region which results from 
analysis of elliptic flow data governed in terrestrial HIC 
experiments 14]. Although the introduced scheme would 
not change, it seems useful to us to repeat these calcu- 
lations under this point of view. It would be interest- 
ing too, to examine the agreement of experimental flow 
data with numerically calculated values explicitly apply- 
ing the KVOR and DD-F models that have passed above 
constraints. 
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We have used here models which do not allow for phase 
transitions. Any possible phase transition that may ap- 
pear in the NSs interior results in a decrease of the max- 
imum NS mass. All the models may well pass the con- 
straint M max > 1.6 M Q (2a uncertainty for PSR J0737- 
3039) but KVR, and even DD-F and KVOR which suc- 
cessfully have passed most of our tests might get prob- 
lems with the restriction M max > 1.9 M@ (la level for 
PSR J0737-3039 and lower limit for 4U 1636-536) if the 
phase transition is sufficiently strong. 

If the phase transition would occur in SNM it would 
also soften the EoS thus modifying the flow constraint 
depicted in Fig. 4. Nevertheless, if the energy gain due 
to the phase transition is not too large the band of the 
flow constraint ^ij is rather broad and all the models 
which already are situated within this band may still 
remain there. However, if we considered the common 
"fiow+maximum mass" constraint as indicated by Fig. 6, 
the NLp and NLpS models could get a problem crossing 
the lower boundary of the thus obtained new band, again 
if the phase transition would be sufficiently strong. 

The charged pion, kaon and the hyperonization tran- 
sition in NSM change the proton and electron concen- 
trations thus affecting the DU threshold. This threshold 
is then pushed up to higher densities. Simultaneously, 
these transitions open new reaction channels, allowing 
for DU reactions which involve new particles: 7r~, K~, 
hyperons and quarks. The appearance of the conden- 
sates and /or filling of the new Fermi seas also affects 
the values of the pairing gaps which are not well known. 
With small gaps the new reaction channels lead to very 
rapid cooling of NS raising the problem to appropriately 
describe the NS cooling. However, if gaps are large these 
reactions might be non-operative and the DU constraint 
may become softer or even non-effective. The threshold 
densities for hyperonization strongly depend on poorly 



known baryon-baryon interactions. In case of repulsion 
the threshold density is pushed up [84| and the situation 
becomes more cumbersome. We avoided studying such 
models here. 

We postpone the analysis of EsoS allowing for phase 
transitions to future work. Besides hyperonization, the 
possibility of a quark matter phase transition should be 
studied. This will be important for both NSs and for the 
planned HICs in the planned CBM experiment at FAIR 
where large baryon densities are to be created. 
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